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Abstract 

Electromagnetic properties of the octet mesons as well as the octet and decuplet baryons are 
augmented in quenched and partially quenched chiral perturbation theory to include 0{a) correc- 
tions due to lattice discretization. We present the results for the SU{3) flavor group in the isospin 
limit as well as the results for SU{2) flavor with non-degenerate quarks. These corrections will 
be useful for extrapolation of lattice calculations using Wilson valence and sea quarks, as well as 
calculations using Wilson sea quarks and Ginsparg- Wilson valence quarks. 
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I. INTRODUCTION 



Lattice gauge theory can provide first principle calculations in the strongly interacting 
regime of QCD, where quarks and gluons are bound into color-neutral hadronic states. These 
calculations, however, are severely limited by the available computing power, necessitating 
the use of light quark masses rUg that are much larger than those in reality. Hence, one 
needs to extrapolate from the quark masses used on the lattice to those of nature. 

A model independent way to do this extrapolation is to study QCD at hadronic scales 
through its low-energy effective theory, chiral perturbation theory (xPT). Since xPT pro- 
vides a systematic expansion in terms of m^/A^, where A^ is the chiral symmetry breaking 
scale, one can, in principle, understand how QCD observables behave as functions of the 
quark mass. In order to address the quenched and partially quenched approximations em- 
ployed by lattice calculations, yPT has been extended to quenched chiral perturbation 
theory (QyPT) B IM IE I E 1^ and partially quenched chiral perturbation theory 




(PQX-PT) li\ 

Recently, we considered the electromagnetic properties of the octet mesons and both 
the octet and decuplet baryons in QxPT and PQxPT [13, EM IE]- Owing in part to the 
charge neutrality of singlet fields, the quenched results are not more singular in the chiral 
limit than their unquenched counterparts. We showed, however, that despite this similarity, 
the quenched results contain singlet contributions that have no analog in xPT. Moreover, 
quenching closed quark loops alters the contribution from chiral logs. For the decuplet 
baryon form factors, for example, quenching completely removes these chiral logs. Many 
others have also observed that the behavior of meson loops near the chiral limit is misrep- 
resented in QxPT, see for example jlE, 21, 2^, 2^, 2^1- On the other hand, PQxPT results 
are devoid of such complications and allow for a smooth limit to QCD. 

Not only are lattice calculations limited to unphysically large quark masses, they are 
also severely restricted by two further parameters: the size L of the lattice, that is not 
considerably larger than the system under investigation; and the lattice spacing a, that is 
not considerably smaller than the relevant hadronic distance scale. To address the issue of 
finite lattice spacing, xPT has been extended (following the earlier work of ^]) in the 
meson sector to 0{a) for the Wilson action j23] and for mixed actions 2^1. Corrections at 
O(a^) have also been pursued 2^, 33|- Corrections to baryon observables have only recently 
been investigated jsj . To consider finite lattice spacing corrections, one must formulate the 
underlying lattice theory and match the new operators that appear onto those in the chiral 
effective theory. This can be done by utilizing a dual expansion in quark mass and lattice 
spacing. Following ji^, |3l[ , we assume a hierarchy of energy scales 



mo -C Ay <C 



1 



(1) 



and ignore finite volume effects. The small dimensionless expansion parameters are 



a A^ 



(2) 



where p is an external momentum. Thus we have a systematic way to calculate 0{a) 
corrections in xPT for the observables of interest. 
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In this work we investigate the 0{a) corrections to the electromagnetic properties of the 
meson and baryon octets, the baryon decuplet, and the decuplet to octet electromagnetic 
transitions in QxPT and PQxPT. We work up to next-to-leading order in the chiral ex- 
pansion and to leading order in the heavy baryon expansion. The paper is structured as 
follows. First, in Section IHl we review PQxPT at finite lattice spacing with mixed actions. 
Since the setup for QxPT parallels that of PQxPT, we will only highlight differences where 
appropriate. Next in Section UTTl we calculate finite lattice spacing corrections to the charge 
radii of the octet mesons to O(e^). This is followed by the calculation of such corrections 
to: the charge radii and magnetic moments of the octet baryons; the charge radii, magnetic 
moments, and electric quadrupole moments of the decuplet baryons; and the decuplet to 
octet electromagnetic transition moments f Sections IIVHVI|) . Corresponding results for the 
above electromagnetic observables in the SU(2) flavor group are presented in Appendix 1X1 
In Appendix ini we determine the 0{a) corrections in an alternative power counting scheme 
for coarser lattices where e ~ cA^. A conclusion appears in Section IVlll 



II. PQxPT AT FINITE LATTICE SPACING 



In partially quenched QCD (PQQCD) H i. Hi [III Q H El [13 the quark part of the 
Symanzik Lagrangian (3^.l33l| to 0{a) is written as 



C = Q{i]p-mQ)Q + aQ a^'^G^, cq Q, 



(3) 



where the second term, the Pauli-term, breaks chiral symmetry in the same way as the quark 
mass term. Here, the nine quarks of PQQCD are in the fundamental representation of the 

and appear in the vector 



graded group SU{Q\3) 

Q = {u, d, s, j, I, r, u, d, s) 
that obeys the graded equal-time commutation relation 



Qnx)Qf (y) - i-ir"'^Qf{y)Qt{^) = 5'^H,,5\-k - y) 



(4) 



(5) 



where a and (3 are spin, and i and j are flavor indices. The remaining graded equal-time 
commutation relations can be written analogously. The different statistics for fermionic and 
bosonic quarks are incorporated in the grading factor 



1 for = 1,2,3,4,5,6 
for /c = 7, 8, 9 



(6) 



The quark mass matrix is given by 

mq = diag(m^., rud, m^, mj, mi, rrir, rriu, rrid, rUs), (7) 
while the Sheikholeslami-Wohlert (SW) [3^ coefficient matrix for mixed actions reads 

CQ = diag(c^ c^ c^ c^ c^ c^ c^ c^ c^. (s) 



If the quark Qi is a Wilson fermion |38|], then {cQ)i = Csw Alternately, if Qi is of the 
Ginsparg- Wilson variety [3^ (e.g., Kaplan fermions ji^] or overlap fermions |41j]), then 
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{cQ)i = 0. Since one expects simulations to be performed with valence quarks that are all 
of the same species as well as sea quarks all of the same species, we have labeled the SW 
coefficients in Eq. (jH)) by valence (v) and sea (s) instead of flavor. In the limit rrij = niu, 
mi = rrifi, and = m<, one recovers QCD at 0{a). 

The light quark electric charge matrix Q is not uniquely defined in PQQCD i^l- By 



imposing the charge matrix Q to have vanishing supertrace, no new operators involving the 
singlet component are introduced. This can be accomplished by ji^ 

Q = diag 1h Qr, Qj, qi, Qr^ ■ (9) 

In addition to the SW term in Eq. Q, the vector-current operator of PQQCD also 
receives corrections at 0{a). There are three operator structures to consider 



0^ = aQQcomQ-f^Q 
0^ = aQQci (iD^'^Q 

0!^ = aD,{QQc2cr^'''Q), (10) 

< — — > 

where = — D^^ and D'^ is the gauge covariant derivative. The form of the matrices 
Co, ci, and C2 in PQQCD is 

c, = diag (cj,cj,cj, 4 44cJ,cJ,c-), (11) 

where and Cj are the coefficients of the vector-current correction operator Oj for valence 
and sea quarks respectively. If the vector-current operator is 0{a) improved in the valence 
(sea) sector, then = ( = ). The operator Oq, which corresponds to a renormalization 
of the vector current, contains a factor of aruQ that renders it 0{e^). Thus contributions to 
electromagnetic observables from Oq are neglected below. The equations of motion which 
follow from Eq. Q can be used to show that the operator O2 is redundant up to 0{a^) 
corrections. Therefore, we need not consider For ease we define the matrix product 
ci,s = Qci. 

A. Mesons 

For massless quarks at zero lattice spacing, the Lagrangian in Eq. (jS)) exhibits a graded 
symmetry SU{6\3)l^SU{6\3)r^U{1)v that is assumed to be spontaneously broken down to 
SU{6\3)v ® U{l)v- The low-energy effective theory of PQQCD that results from expanding 
about the physical vacuum state is PQxPT. The emerging 80 pseudo-Goldstone mesons can 
be described at O(e^) by a Lagrangian which accounts now for the two sources of explicit 
chiral symmetry breaking 



C = y str (D^StD^S) str (mgE + m^E^) +aA„ str (cqE + c^S^) +a^f'^o^^,^o-^4,^l 

(12) 

where 

E = exp(^)=5^ (13) 
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^-{"H). (14) 



X M 

f = 132 MeV, and the gauge-covariant derivative is -D^S = + ieAfj,[Q, S]. The str() 
denotes a graded flavor trace. The M, M, and x are matrices of pseudo-Goldstone bosons 
and pseudo-Goldstone fermions, see, for example, [i^. Expanding the Lagrangian in ()12|) 
one finds that to lowest order mesons with quark content QQ' have mass 



„2 
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^QQ' = -p [^m{.mQ + mg/) + aA„(cQ + Cq,)] . (15) 



The flavor singlet field is $o = str($)/v^- It is rendered heavy by the U{\)a anomaly 
and can be integrated out in PQxPT, with its mass /iq taken on the order of the chiral 
symmetry breaking scale, /xq — ^ ^x- -'■^ ^^^^ limit the propagator of the flavor singlet field is 
independent of the coupling a and deviates from a simple pole form ^ 3|. In QxPT, the 
singlet must be retained. 



B. Baryons 

In PQQCD there are baryons with quark composition QiQjQk that can contain all three 
types of quarks. The spin-1/2 baryons are embedded in the 240-dimensional super-multiplet 
Bijk, that contains the familiar octet baryons, while the spin-3/2 baryons are embedded in 
the 138-dimensional super-multiplet that contains the familiar decuplet baryons [6t. i4d| . 

At leading order in the heavy baryon expansion and at 0{a), the free Lagrangian for the 
Bijk and T^jj^ fields is given by ja |31 



C = i{Bv VB) + 2aM (BBM+) + 2/3^ (BM+B) + 2aM (BB) str {M+) 
+2aA {BBA+) + 2Pa (BA+B) + 2(Ta (BB) str {A+) 
-I (r% ■ VT,) + A (T%) + 27M (T^M+T,) - 2a u {T%) str 
+2jA (J^A+T;) - 2a A {r%) str {A+) , (16) 

where = | {S^^nqS} + ^mqC,) and A+ = |a (C^cq^^ + ^cq^). Here A is the mass split- 
ting between the 240 and 138 in the chiral limit. The parenthesis notation used in Eq. ()16|) 
is defined in so that the contraction of flavor indices maintains proper transformations 
under chiral rotations. Notice that the presence of the chiral symmetry breaking SW op- 
erator in Eq. (jS)) has lead to new 0{a) operators (and new constants a a, Pa, <^a-, 1a-, and 
'a a) in Eq- (fT^ - The Lagrangian describing the interactions of the Sjjfc and T^-^ with the 
pseudo-Goldstone mesons is 



L = 2a {BS^BA;)^2P {BS^A^B)+2n {T"S^A^%) + ^C [{T" A,B) + [BA.T")] . 



(17) 



The axial- vector and vector meson fields A^^ and are defined by: A^^ = ^ {C.d^C,'' — ^^d'^^) 
and ¥^" = 1 {^d"^^ + ^^df"^). The latter appears in Eq. (fTBj) for the covariant derivatives of 
Bijk and Tijk that both have the form 

{V'^B)ijk = d^Bi,k + {V>')iiBi,k + H'''^''^'''''^\V>'),^B,rnk + {-Y''^^^^^^ (18) 
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The vector 5''^ is the covariant spin operator 4^ 4^ i^l- The parameters C and 7i entering 
in Eq. ()17|) are identical to those in QCD, while the parameters a and (3 can be related to 
the familiar constants D and F of QCD 

a = -D + 2F and (3 = --D + F. (19) 
3 3 

In QQCD these identifications cannot be made. 

III. OCTET MESON ELECTROMAGNETIC PROPERTIES 

The electromagnetic form factor G{q^) of an octet meson has the form 

mp')\j^mp)) = G{q'){p+pr (20) 

where q'^ = {p' — p)^. At zero momentum transfer G{0) = Q, where Q is the charge of 0. 
The charge radius r is related to the slope of G(g^) at = 0, namely 

< >= 6 • (21) 

Recall, at one-loop order in the chiral expansion the charge radii are O(e^). 
There are two finite-a terms in the 0{e^) Lagrangian [29| 



C = aAA^MD.T."^ D^Y.)McQ^ + c^^t) + a^,5^str(Z}^StZ}'^S(cQS + c^St)) (22) 



that contribute to meson form factors at tree level. The new parameters aA,i and a^.s in 
Eq. (j22p are finite lattice spacing analogues of the dimensionless Gasser-Leutwyler coeffi- 
cients 0:4 and as of xPT |4S|. The above terms contribute to meson form factors at 0{e^) but 
their contributions are independent of q^ and annihilated by the corresponding wavefunction 
renormalization, thus ensuring charge non-renormalization. 

The SW term can potentially contribute at 0{t^) when A+ is inserted into the kinetic 
term of the leading-order C in Eq. (fT^ . Contributions to form factors from such terms vanish 
by charge non-renormalization. Insertions of Aj^ into the ctg term of the Gasser-Leutwyler 
Lagrangian produces the 0{e^) terms 



C = 



imiA^F^, str {{Q+, A+jD^'^D'^^ + {Q+, A+jD^"^^ D'^) 
+im2A^F^u str {Q+D^'^A+D''^^ + Q+D^SU+D'^S) 
+im3A^F^^str (Q+D'^SD'^St + Q+D'^S'^D'^S) str(^. ^ 



(23) 



where we have defined Q+ = | (^^Q^^ + ^Q^). These terms contribute at 0{e'^) to the 
charge radii and can be ignored (see Appendix ^ for discussion relating to larger lattice 
spacings) . 

Additionally we must consider the contribution from the vector-current correction oper- 
ator Oi in Eq. (ITUI). In the meson sector, the leading operators in the effective theory 
can be ascertained by inserting aA^ci^g in place of Q in the operators that contribute to 
form factors. The effective field theory operators must also preserve the charge of the meson 
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0. It is easiest to embed the operators in a Lagrangian so that electromagnetic gauge 
invariance is manifest. To leading order, the contribution from Oi is contained in the term 

C = iaA,9 aA^F^^ str {ci^Qd^Wi:^ + ci^Qd^'E^d^i:) . (24) 

Thus the correction to meson form factors from Oi is at 0{e^). 
The charge radius of the meson to C(e^) then reads 

< >= + — i— Ax log (25) 

where X corresponds to loop mesons having mass mx [the masses implicitly include the 
finite lattice spacing corrections given in Eq. (jl5|) . otherwise the expression is identical to 
the a = result]. The coefficients Ax in PQxPT appear in Ref. In the case of QxPT, 
the coefficients Ax = for all loop mesons and there are no additional contributions from 
the singlet field at this order. Thus there is neither quark mass dependence nor lattice 
spacing dependence in the quenched meson charge radii at this order. 



IV. OCTET BARYON ELECTROMAGNETIC PROPERTIES 

Baryon matrix elements of the electromagnetic current J'^ can be parametrized in terms 
of the Dirac and Pauli form factors Fi and F2, respectively, as 

(B{p') \J^\ B{p)) = u{p') i^v^F,{q') + '^^^q^F,{q')^u[p) (26) 

with q = p' — p and Mb is the degenerate octet baryon mass. The Dirac form factor is 
normalized at zero momentum transfer to the baryon charge: -Fi(O) = Q. The electric 
charge radius < r\>, and magnetic moment ^ can be defined in terms of these form factors 
by 

and 

/i = F2(0). (28) 

Recall, that the one-loop contributions in the chiral expansion to the charge radii are 0{e^), 
while those to the magnetic moments are 0{e). 

There are no finite-a operators in Eq. that contribute to octet baryon form factors. 
As in the meson sector, however, the SW term could contribute when A+ is inserted into 
the Lagrangian. Here and henceforth we do not consider these insertions into the kinetic 
terms in Eq. (fT^ because their contributions alter the baryon charges and will be canceled 
by the appropriate wavefunction renormalization. 

The SW term, however, does contribute when Aj^ is inserted into the charge radius and 
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magnetic moment terms. For the charge radius, we then have the 0{a) terms 

Ay- 

X 



1 

a; 

h 

A. 



^x 

+_L {BBQ+) + h (BQ+B)] v^d.F^^' sti{A+) 
^x 

+^ {BB) v.d^F^" str(Q+^+), (29) 

that contribute at 0{e^) to the charge radii and are thus neglected. Insertions of A+ into 
the magnetic moment terms produce 



C = ih\{-)^''^+''^'^^'''^+''^'^W''[S^,S,]{Q+,A+Y^'B''^' F^"^ 
+ih'^ W'\S^, S,]{Q+, A+Y'B''^ F^" 
+ih'^{-f^'^'^^+'^^'^W^'[S^, SJ\ Q^A^'B''^'''F>"' 

+t [fo'e {B[S„S,]BQ^)+b', {B[S„S,]Q+B)] F^''sti{A+) 

+1 fe's {B[S^, S,]B) F^" str(Q+ A), (30) 



which are 0{e^) corrections to the magnetic moments and can be discarded |31|. 

Finally we assess the contribution from the operator in Eq. (|lUj). As in the meson 
sector, the charge preserving operators can be constructed by the replacement Q —>■ aA^ci^Q 
in leading-order terms. Again it is easier to embed these operators in £ so that gauge invari- 
ance is transparent. For the charge radius, the leading vector- current correction operator is 
contained in the term 

[cA,a {BBc,,q) + CA,I3 {Bc,,qB)] V^d^F'^'', (31) 

which leads to 0{e^) corrections. For the magnetic moment operator, such a replacement 
leads to 

f} ft 

£ = - [fiA,a {B[S„ S,]Bc^,q) + fiA,p {B[S^, S,]c,,qB)] F^\ (32) 

and corrections that are of higher order than the one- loop results See Appendix iBl for 
results in an alternate power counting scheme. 
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To 0{e'^) the baryon charge radii are thus 



< r| > 



^ ((5c_ + aDC+) + T^iQl^F + ttD/iZ)) 



2M' 



B 



X 



Ax log 



m 



X 



5 /?x log 



m 



X 



10/?^e?(mx,A,/i) 



(33) 



and the magnetic moments to 0{e) read 



fi = {Q ^ip + aD /Ud) + ^ + Px^{mx, A, /i)] 



47r/ 



(34) 



The a-dependence is treated as implicit in the meson masses. The PQxPT coefficients Ax, 
Px, and P'x can be found in [13, 13 along with the functions J-'{mx, A, /i) and Q{mx, A, /i). 
The quenched charge radii at 0{e^) are similar in form (although Ax = in QxPT) due 
to the lack of singlet contributions at this order. The QxPT coefficients and appear 
in 2^1. The quenched magnetic moments, however, receive additional contributions from 
singlet loops. The relevant formula of 2^ are not duplicated here in the interests of space 
but only need trivial modification by taking into account the a-dependence of meson masses. 



V. DECUPLET BARYON ELECTROMAGNETIC PROPERTIES 

Decuplet matrix elements of the electromagnetic current J'' can be parametrized as [ii 

{T{p')\r\T{p)) = -u,{p')0^^''u,{p), (35) 

where w^(p) is a Rarita-Schwinger spinor for an on-shell heavy baryon. The tensor O'^'"^ can 
be parametrized in terms of four independent, Lorentz invariant form factors 



where the momentum transfer q = p' — p. The form factor Fi{q^) is normalized to the 
decuplet charge in units of e such that -Fi(O) = Q. 

The conversion from the covariant vertex functions used above to multipole form factors 
for spin-3/2 particles is explicated in [i^. For our calculations, the charge radius 

< r| >= 6 - TTT^ [2g - 3F2(0) - Gi(0)]l , (37) 



the magnetic moment 
and the electric quadrupole moment 



dg2 12M^ J 

/x = F2(0), (38) 



\Gm- (39) 



The charge radii are 0{e^) at next-to- leading order in the chiral expansion, while the mag- 
netic moments are 0{e) and the electric quadrupole moments are 0{e'^). At one-loop order 
in the chiral expansion, the magnetic octupole moment is zero. 
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There are no finite-a operators in Eq. p6p that contribute to decuplet baryon form factors. 
The SW term can potentially contribute when is inserted into the Lagrangian. There 
are three such terms: the charge radius, magnetic moment, and electric quadrupole terms. 
Insertions of A+ into the charge radius term produces 

+ ^ (T'^Q+X) v^d^F'^" str(A) + ^ (T'^X) v^d.F^" str(Q+^+). (40) 

These contribute to decuplet charge radii at 0{e^). As in the octet sector, insertions of A+ 
into the magnetic moment term, namely 

+1 4 (r^Q+X) F^" str(^+) + 1 4 {T^%) F^" str(Q+ A), (41) 

produce 0{e^) corrections. Likewise, insertions of A+ into the electric quadrupole term have 
the form^ 

C = ^T^^'^^'{Q+,A+Y'T''^''''^v''d^F,^ + ^{-fi'^'^^+'^^'^^ 
^x ^x 

+ §- {T^'^Q+T''^) v'^d.F,^ str(^+) + ^ (T^/^T'^}) v^d^F,^ str(Q+A), (42) 
^^x ^^x 

and produce O(e^) corrections. All of these corrections are of higher order than the one- loop 
results. 

Finally we assess the contribution from the operator Oi in Eq. (fTUI) . The effective op- 
erators can be constructed by replacing Q by aA^ci^g in leading-order terms. Embedding 
these terms in a Lagrangian, we have 

C = ^ {r''c,,Q%) v.d^F^'^ + 3^a/.^,, [T,c,,Q%) F^^ - [T^^c,,Qr^) v'^d.F^^. 

X X 

(43) 

Each of these terms leads to corrections of higher order than the one-loop results and can 
be dropped. Thus at this order the only finite lattice spacing corrections to decuplet elec- 
tromagnetic properties appear in the meson masses. For reference, the expressions are 



5^ Ax log -f - - -— ^ ^xQimx, A, /i), (44) 



3 16n^P ^ " ° /i2 27 Wtt^P ^ 



/i = 2/ic(5 - gg^^2 Ax^("^x, A, /i) - g^^^ 2^ v4xmx, (45) 



^ The action of '^ '^ on Lorentz indices produces the symmetric traceless part of the tensor, viz.. 
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and 



(46) 

The coefficients Ax are tabulated in [3]. Extending the resuh to QxPT, where Ax = 0, 
one must include additional contributions from singlet loops. With finite lattice spacing 
corrections, the expressions are identical to those in [l^ except with masses given by Eq. (|TK|) . 
Thus for brevity we do not reproduce them here. 



VI. DECUPLET TO OCTET ELECTROMAGNETIC TRANSITIONS 

The decuplet to octet matrix elements of the electromagnetic current appear as [Ji 

{B{p)\r\T{p')) = u{p)0^%{p'), (47) 

where the tensor O^^ can be parametrized in terms of three independent, Lorentz invariant 
form factors 

O^^ = {q . Sg^^ - SV) ^ +{q-vg^^- v^q^) q ■ S^^fl 

' Mb ^ [2MbY 

+ {<f9''-<l'q')S-q^l (48) 

where the photon momentum q = p' —p. At next-to-leading order in the chiral expansion, we 
recall that Gi{q^) is 0(e) while G2{q^) and G^{q^) are 0(e°). The conversion of these vertex 
covariant functions to multipole form factors is detailed in [sOl • The multipole moments up 
to C(e) are^ 

- es;) + (5 + (4) "^-^'^ ^ l"^^"^- 

There are no new finite-a operators in Eq. ()16|) that contribute to decuplet to octet tran- 
sition form factors. Insertion of A+ into leading-order transition terms leads to corrections 
of O(e^) or smaller. For completeness the terms are: 



+ih F^" str(^+), (52) 



Here, we count e ~ A/Mb- 
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for the magnetic dipole transition; and 

^x 

^x 

+^ (BS^^Q+r''^) v^d^F^o. str(^+), (53) 

for the quadrupole transition. Finally, insertion of A+ into the PQxPT term proportional 
to i(BS''Q.T'')d'^d^F^a leads to 

C = '^^B^^'S^QlA%T^^^d''d>'F\+'^^B''^'S^J^^ 

X x 

X 



A2 

X 



+^ (i35^Q+r.) a"a'^F'^„str(^+), (54) 



A2 

X 



for the Coulomb quadrupole transition. 

Similarly, constructing in the effective theory by replacing Q with aA^ci^Q in the 
transition operators leads to terms of at least 0{e^) which are contained in the terms 



JO- - ta^iA,T^{BS,Cr,Q%)F^'' + ^^^^{BS^^Cr,QT''^)v"d,^^^ 

'■^{BS>^c,,QT^)d"d,F,^. (55) 



A2 

X 



All of these corrections from effective Oi operators are of higher order than the one-loop 
results. Thus at this order, the only finite lattice spacing corrections to the transition 
moments appear in the meson masses. For reference the expressions are 



+47rC(Z^-F)-f / dx{l-x)J^{mx,-xA,i^), 

\ X -^0 



(56) 



-16C{D - F)J2(^x dxx{l-x)g{mx,-xA,iJ,)^, (57) 
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and 



^3(0) = -16-r|-^ / dxx{\-x)[x 




\\ Amx 
2) m\ - x2A2 



+ CiD-F)(3^'JZ 




(58) 



The coefficients Pb and Pt are tabulated in [l9| along with the function lZ{x). Extending 
the result to QxPT, where the coefficients are replaced with their quenched counterparts 
and /5y , one must include additional contributions from singlet loops. With finite lattice 
spacing corrections the expressions are identical to those in [19] except with masses given 
by Eq. (|T5|l . Thus for brevity we do not reproduce them here. 

VII. CONCLUSIONS 

Above we have calculated the finite lattice spacing corrections to hadronic electromagnetic 
observables in both QxPT and PQxPT for the SU{'i) flavor group in the isospin limit and 
the SU{2) group with non-degenerate quarks. In the power counting scheme of j29[ l3ll |. 
0{a) corrections contribute to electromagnetic observables at higher order than the one- 
loop chiral corrections. Thus finite lattice spacing manifests itself only in the meson masses 
at this order. 

In practice one should not adhere rigidly to a particular power-counting scheme. Each 
observable should be treated on a case by case basis. The actual size of a and additionally the 
size of counterterms are needed to address the relevance of 0{a) corrections for real lattice 
data. For this reason we have presented an exhaustive list of 0{a) operators relevant for 
hadronic electromagnetic properties. In an alternate power counting for a coarser lattice (as 
explained in Appendix^, some of the operators listed above contribute at the same order 
as the one-loop results in the chiral expansion. The corrections detailed in Appendix FBI in 
the baryon sector are also necessary if one goes beyond the heavy baryon limit and includes 
recoil terms to the Lagrangian, even in the original power counting of Eq. (j21). 

Knowledge of the low-energy behavior of PQQCD at finite lattice spacing is crucial to 
extrapolate lattice calculations from the quark masses used on a finite lattice to the physical 
world. The formal behavior of the PQQCD electromagnetic observables in the chiral limit 
has the same form as in QCD. Moreover, there is a well-defined connection to QCD and one 
can reliably extrapolate lattice results down to the quark masses of reality. For simulations 
using unimproved lattice actions (with Wilson quarks or mixed quarks), our results will aid in 
the continuum extrapolation and will help lattice simulations make contact with real-world 
data. 
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APPENDIX A: HADRONIC ELECTROMAGNETIC PROPERTIES IN SU{2) 
FLAVOR WITH NON-DEGENERATE QUARKS AT FINITE LATTICE SPACING 



In this Appendix, we consider the case of SU{2) flavor PQQCD^ and summarize the 
changes needed to determine finite lattice spacing corrections to the electromagnetic prop- 
erties of hadrons considered above. For the two flavor case, we keep the up and down valence 
quark masses non-degenerate and similarly for the sea-quarks. Thus the quark mass matrix 
reads 

mg^^^^ = diag(m„, rud, mj, mi, ruu, rud), (Al) 

while the SW matrix is 

= diag(c^ c^c^c^c^c^). (A2) 

Defining ghost and sea quark charges is constrained only by the restriction that QCD 
be recovered in the limit of appropriately degenerate quark masses. Thus the most general 
form of the charge matrix is 

= diag 0, g„ qi, g„ , (A3) 

which is not supertraceless. Analogous to the three flavor case, the vector-current will receive 
0{a) corrections from the operators in Eq. (fTUI) of which only the operator Oi is relevant. 
The coefficient matrix associated with this operator is 

cf^''> = dis.g{cl,cl,clcl,clc\). (A4) 

The 0{a) operators listed above in Sections IIIIUVII are the same for the SU{2) flavor 
group, however, the coefficients have different numerical values. Additionally there are 
operators involving str(Q^'^*^^''). These are listed for each electromagnetic observable below. 

Octet mesons 

In the meson sector, one has the additional term 

£ = im^A^F^^sti {A+D^'^D"^^ + A+D^'^^D"^) str(Q+^^^^). (A5) 

Octet baryons 

In the octet baryon sector, there are terms which originate from A+ insertions 
£ = [69 {BBA^)+ho (BA+B)] v^d^F'^'' sti{Qf^'^) 

+^ {BB) v^d,F^''stT{Qf"''^)stT{A+) 
^x 

+t [b', {B[S^, S,]BA+) + b[, {B[S^, SM+S)] str(Qf (')) 

+ih\, {B[S,, S,]B) F'^'^str(Qf (')) str(^+), (A6) 



For brevity we refer to SU{A\2) PQQCD as SU{2). The distinction will always be clear. 
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and additional vector-current correction operators 

(A7) 



C = ^ v.d^F'^'' str(Q^^(2)cf^'0 + (^[^m- S^^) i^'^^str(Q^^(2)^6i 



Decuplet baryons 

Next in the decuplet sector there are terms that result from insertions 

C = ^(rM+T^)t;^9,F^'^str(Qf(')) + ^(r-r.)t;^9,F^^str(Qf('))str(^+) 
^^x ^^x 

+z 4 F^-" str(Qf + 1 4 {T,%) F^'^ str(Qf str(^+) 

+ ^ (r^M+T'^}) str(Qf + ^ (T^'^T'^}) str(Qf str(^+) 

^^x ^^x 

(A8) 

and also further vector-current correction operators 

£ = (T'^T.) t;,9.F'^'^str(Q^^(2)cf (^)) + 3m /i'^,^ (T,T.) F'^'^ str(Q^^(2)cf (^)) 



X 



.^^^ (^{.^^1) ^a^^^^^ str(Q^^(2) cf(')). (A9) 



Ax 



Baryon transitions 

Finally for the transitions, there are only new A+ insertions 

£ = tte (BS,A+%) F^'^str(Qf + ^ (BS^'^A+T''^) v^d^F,^stT{Qf'^'^) 

^^x 

+^ (BS^A+%) d^d^F\ str(Qf (AlO) 

X 

For each electromagnetic observable considered above, contributions from all 0{a) op- 
erators in the effective theory are of higher order than the one-loop results in the chiral 
expansion. Thus one need only retain the finite lattice spacing corrections to the meson 
masses and use the previously found expressions for electromagnetic properties in SU{2) 

pQxPT HI El Ei, Sim. 



APPENDIX B: COARSE-LATTICE POWER COUNTING 

In this Appendix, we detail the 0{a) corrections to electromagnetic properties in an 
alternate power-counting scheme. We imagine a sufficiently cursed lattice, where aA^ can 
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be treated as C(e), so that"^ 




■X' 



(Bl) 



In this case, there are known additional O(a^) corrections j29| to the meson masses that are 
now at O(e^) and must be included in expressions for loop diagrams. The free Lagrangian 
for Bijk and 7^^^ fields contains additional terms of C(a^) that correct the baryon masses, 
and modify the kinetic terms. Potential contributions due to the latter, whatever their form, 
must be canceled by wavefunction renormalization diagrams. The only contribution of 0{a^) 
could come from tree-level electromagnetic terms but these are necessarily higher order. 
Thus in this power counting there are no unknown O(a^) corrections for electromagnetic 
properties. 

The only possible corrections come from the 0{a) operators assembled above. A few of 
these do contribute at tree level and are spelled out below. 

Octet mesons 

The 0{a) corrections to the meson form factors are now 0{e^) in the power counting. 
While the meson charge radii at NLO in the chiral expansion are at 0{e^), further corrections 
in the chiral expansion are at 0{e'^). Thus one can use the 0{a) operators to completely 
deduce the charge radii to 0{e^) [apart from 0{e^) corrections to the meson masses]. These 
0{a) operators are given in Eqs. (j^Hj) and ()24|) and yield a correction S < r'^ > to the meson 
charge radii of the form 



Notice that there are no corrections associated with an unimproved current operator in the 
sea sector since is absent. 

In the case of SU (2) flavor, there is an additional contribution from the operator in 
Eq. ()A5|) . At tree level, however, this operator vanishes. The only correction to Eq. ()B2|) in 
changing to SU (2) flavor is to replace 3c* with 2c'^ which reflects the change in the number 
of sea quarks. 

Octet baryons 

For the octet baryon electromagnetic properties, the 0{a) corrections to the charge radii 
are now 0{e^) and can be dropped as they are the same order as neglected l/Mg corrections. 
The magnetic moments, however, do receive corrections from local operators. Specifically, 
the 0{a) operators which contribute to magnetic moments at 0{e) are insertion of into 
the magnetic moment operator given in Eq. (j3U|) and corrections given in Eq. (|32j). 



^ This power counting coupled with the chiral expansion is most efficient for valence Ginsparg- Wilson quarks 
where 0{a) corrections vanish. We thank Gautam Rupak for pointing this out. 



5 <rl>=Q 



24aA. 



■X 



[d'{2mi + + 3c''m3 + claA,9] 



(B2) 



16 



TABLE I: The coefficients 



A and B for the octet baryons. 



P 
n 

S+ 
SO 
A 

S^A 

s- 



A 


B 


1 





1 
i 


1 
i 


3 


3 


1 





1 


1 


6 

1 


6 

1 


6 


6 


1 


1 


2\/3 


2V3 


2 


1 


3 


3 


1 


1 


3 


3 


2 


1 


3 


3 



Calculation of these corrections yields a shift in the magnetic moments 



dfi = aMs I A (^b[ + h'^^ - B (26'2 + - b',) + 3c' b'^ - B b'^^ 



+ Cic'-cnq,irb's + ^ 



(B3) 



where qjir = qj + qi + qr- The coefficients A and B are listed for octet baryons in Table 
m while C = 1 for all octet magnetic moments and C = for the AE° transition moment. 
Notice that there are no corrections associated with an unimproved current operator in the 
sea sector. 

In the case of SU{2) flavor, there are additional contributions given in Eqs. flA6|) and 
rfj). For the proton and neutron, we have 



A{b[ + -6; - B i2b', + 6'3 - b',) + -ib', + b[,) 



+ 2c'{^Ab',^Bb', + h[, 



+ 



(B4) 



where qji = qj + qi. 



Decuplet baryons 

For the decuplet baryon electromagnetic properties in coarse-lattice power counting, the 
0{a) corrections to the charge radii are C(e^) and the corrections to the electric quadrupole 
moments are 0(e), both of which are higher order than the one-loop results. The mag- 
netic moments, however, do receive corrections from local operators. Specifically, the 0{a) 
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operators which contribute to magnetic moments at 0{e) are A+ insertions into the mag- 
netic moment operator given in Eq. (PT|) and Of correction operators given in Eq. ()43p. 
Calculation of these corrections yields a shift in the magnetic moments 



Sfi = 2aMB 



-c'^Q (24 + 4) + c'Q 4 + (c^ - c")g,7.4 + clQ ^a. 



(B5) 



Notice that in 5't/(3) strQ = 0, hence there is no dependence on c\ in the above result. 

In the case of SU{2) flavor, there are additional contributions given in Eqs. dMl) and 
(IA9jl . The corrections to the A quartet magnetic moments are then 



^ 2aMB { ^c" (2Q d[ + Q 4 + 4) + ^c' {Q 4 + 4) + 



d' 



+ cl [Q fiA,c + (1 - 3gj7)/i!4,7] + 



(B6) 



Baryon transitions 

For the decuplet to octet electromagnetic transitions in coarse-lattice power counting, 
the 0{a) corrections to (^2(0) and 6*3(0) are 0{e) which are of higher order than the one- 
loop results. The Gi(g^) form factor does, however, receive corrections from local operators. 
Specifically these 0{a) operators which contribute to Gi(0) at (9(e) are the insertions of 
A+ into the magnetic dipole transition operator given in Eq. ()52|) and the vector-current 
corrections given in Eq. ()55|). Calculation of these corrections yields a shift of Gi(0) 

5^1(0) = aMB ttT^I jc'' (^h + t2 + h- + 3c' t, + cl /^^.T^lj , (B7) 
where the transition coefficients ar appear in Again, at this order the result is inde- 



pendent of 0{a) improvement to the electromagnetic current in the sea sector. In the case 
of SU(2) flavor, there is an additional dipole operator given in Eq. flMOll . At tree level, 
however, this operator vanishes. The only correction to Eq. flB7|) in changing to SU{2) flavor 
is to replace 3c* with 2c* which reflects the change in the number of sea quarks. 
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